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Abstract 

The matrix model of topological field theory for the moduli space of p-th 
spin curves is extended to the case of the Lie algebra of the orthogonal group. 
We derive a new duality relation for the expectation values of characteristic 
polynomials in the antisymmetric Gaussian matrix model with an external 
matrix source. The intersection numbers for non-orient able surfaces of spin 
curves with k marked points are obtained from the Fourier transform of the 
k-point correlation functions at the critical point where the gap is closing. 
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1 Introduction 



For Riemann surfaces, it is known from Witten's conjectures [2] and Kont- 
sevich's derivation [1] that the intersection numbers of the moduli space of 
curves with marked points may be obtained from an Airy matrix model. Fur- 
thermore, higher Airy matrix models have been shown to give the intersection 
numbers of the moduli space for p-spin curves [3]. 

Recently, a duality relation has been applied to this problem [4, 5, 6]. The 
derivation relied on a duality between the higher p-th Airy matrix models 
and Gaussian matrix models in an external matrix source at critical values 
of this source [7, 8] ; the intersection numbers are then easily obtained from 
this dual model. When p=2, the model reduces to Kontsevich's model ; its 
dual connects to the behavior of correlation functions near the edge of the 
semi-circle spectrum [9, 4]. 

The moduli space of p-th spin curves is described by random Hermitian 
matrices, the Lie algebra of the unitary group U(N). It is of interest to 
extend this moduli space of spin curves for non-orientable surfaces, both in 
the fields of open string theory and of quantum chaos. 

For obtaining non-orientable surfaces from standard loop-expansions of 
matrix models, a first possibility would be to use real symmetric matrices. 
The Euler characteristics of the moduli spaces of real algebraic curves with 
marked points may be obtained from the real symmetric matrix model [10]. 
However, for the intersection numbers of the moduli space of curves with 
marked points, this real symmetric matrix model remains difficult to solve 
when one extends Kontsevich Airy matrix model to non-orientable surfaces 
[11]. 

In this article we have chosen, instead of real symmetric matrices, to 
consider real antisymmetric matrices, the Lie algebra of the SO(N) group 
(we assume that N is an even integer). Let us note that the Gaussian random 
matrix model of the classical groups O(N) and Sp(N) appeared earlier in 
the literature in the studies of the moments of the /^-functions [12] and in 
the study of the spectrum of excitations inside superconducting vortices [13]. 

For the O(N) matrices, there is a HarishChandra formula for the integrals 
over the orthogonal group. Thanks to this integral formula, which is similar 
to the unitary case, generating functions of the intersection numbers for non- 
orientable surfaces become calculable. We shall discuss a duality relation for 
the O(N) case, which is surprizingly similar to the U(N) case ; then we 
compute explicit expansions for the Fourier transforms of the correlation 
functions of the dual models, and obtain the intersection numbers. This 
study may shed some light on the moduli space of curves on non-orientable 
surfaces. 
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2 Duality relation 



Let us first state the basic duality relation which will be used in this article. 
Theorem 1: 



N 



< [] det(A a • I - X) > A =< J] det(a n • I - Y) > A 



(2.1) 



01=1 



n=l 



where X is 2N x 2N real antisymmetric matrix (X* = —X ) and Y is 
2/c x 2/c real antisymmetric matrix ; the eigenvalues of X and Y are thus 
pure imaginary. A is also a 2iV x 2N antisymmetric matrix, and it couples 
to X as an external matrix source. The matrix A is 2k x 2k antisymmetric 
matrix, coupled to Y. We assume, without loss of generality, that A and A 
have the canonical form : 
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(2.2) 



i.e. 



A = a±v 
A is expressed also as 



(InV, v = iv i = 



A = Ai« © • • • © \ k v. 



1 
-1 



The characteristic polynomial det(A-I— X) has the 2N roots (±iAi, 
The Gaussian averages in (2.1) are defined as 



(2.3) 

(2.4) 
,±iA n ). 



< ■ • • >a = — / dXe^ 



±trX 2 +trXA 



< • • • >A 



Z A 
-I 

z a 7 



dYe* 



itrY 2 +trYA 



(2.5) 



in which X is a 2N x 2N real antisymmetric matrix, and Y a 2k x 2k real 
antisymmetric matrix ; the coefficients Z\ and are such that the expec- 
tation values of one is equal to one. The derivation of Theorem 1 relies 
on a representation of the characteristic polynomials in terms of integrals 
over Grassmann variables, as for the U(N) or U(N)/0(N) [12, 14, 15] cases. 
Given the complexity of the intermediate steps of the derivation for the O(N) 
case, the simplicity of the result is striking. The derivation is given in ap- 
pendix A. 
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3 Higher Airy matrix models 



From the theorem 1, we can obtain easily the higher Airy matrix models. We 
first consider the simple case in which the source A is a multiple of identity 
: a„ = a. Then we write 



JV 



< J]det(a n -I-y) > A = < [det(a • I — Y)] N >\ 
i=i 

= _L f ^y e iVtrlog(a-I-y)+itrr 2 +tryA ^.\) 

Z A J 

Expanding the logarithmic term, and noting that the traces of odd powers 
of Y vanish since Y is antisymmetric, we obtain 

< [det(a - Y)] N > A = ±- [ dYe 2km °* a -^-^ Y2 -^ Yi +--- +m (3.2) 

Z a J 

Chosing a 2 = JV, the coefficient of trF 2 vanishes. Then one rescales Y — > 
NiY, and A — > N~*A. After these rescalings, we obtain in the large N limit, 
the higher Airy matrix model, 

Z = J dYe-* tlY4+t * YA (3.3) 

Note that higher powers of Y 2n disappear in this scaling limit since they are 
given by 

■ N 2 -rtrY 2n ~ A-t +1 trF 2ra (3.4) 

which vanish in the large-N limit for n > 2. 

By appropiate tuning of the a n 's, and corresponding rescaling of Y and 
A, one may generate similarly higher models of type 

Z = J dYe-^ trYP+1+trYA (3.5) 

where p is an odd integer. These models are similar to the generalized Kont- 
sevich model in the unitary case, which gives the intersection numbers of the 
moduli spaces of p-th spin curves. However, the matrix Y being real and 
antisymmetric, the partition function Z is very different from the unitary 
case and non-orientable surfaces lead to different intersection numbers. 



4 Expansion in inverse powers of Lambda 

The free energy F = \ogZ can be expanded in powers of trA~ m as in the 
unitary case. This is done through the Harish Chandra formula [16] for the 
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integration over the orthogonal group g = SO(2N). We may take Y and A 
in canonical form (2.1) without loss of generality : then the Harish Chandra 
integral reads [16] 

Theorem 2 (Harish Chandra) : 

N 

E (detw)exp[2 £ w(y j )X j ] 
( dge^ Y ^ = C weW ^ , i=\ 2 x2 . (4.1) 

JSO(2N) n (^-^)(A, 2 -A2) 

l<j<k<N 

where C = (2iV-l)! IIj=f ^j-l)!, and W is the Weyl group, which consists 
here of permutations followed by reflections (z/i — > ±yi ;i = 1, • • • , N) with 
an even number of sign changes. 

In the appendix B, we give a more explicit determinantal expression for 
this HarishChandra integral for the orthogonal group SO(2N). 

For the p = 3 critical model (3.5), let us compute the perturbation ex- 
pansion of the free energy F = logZ. From Th. 2, we obtain (see appendix 
B) 

Z = f A dv n(y?2 ~ ^ e-5 E«?-2E (42) 

To obtain the series in powers of 1/A, we change A« — > Af and make a shift 
Hi — ► 2/i — Aj to eliminate the y^Af term. Then, the problem reduces to cubic 
and quartic perturbations with a Gaussian weight. For instance, k = 2 for 
0(2k) case, we have from Th.2, 



/ 



,2 „,2 



Ai — A 2 

1/1 1 x 1,1 F, ^F / 
1 + 72<Af + Al> + 12<A? + Al> + (4 ' 3) 



where we have dropped a normalization constant. In that simple k = 2-case 
it is easy to make this calculation directly without use of Th.2, and the results 
do agree. The coefficients ^ and ^ are universal factors independent of k. 
For general k, it is useful to write the Vandermonde product as a determinant 
with the yf, and the evaluation is straightforward. We introduce the following 
parameters, similar to those of the unitary case, 

j— p— n(p+2) ^ J" jj 

t nJ = (p)^n^ JJ (l p + j + 1) £ -t^jtt (4.4) 

Z=0 i=l A i 

with here p=3. Note that the normalization of the t n j is slightly different 
from the unitary case, in which p is replaced by — p in the first factor of 
(4.4). However, this definition is not approriate for the half-integer genus 



4 



(non orientable surface), since the first factor becomes irrational number. 
So, we define in this half-integer genus case (condition will appear below in 
(4.9)) as 

n— 1 k i 

^=n('p+j+ i )ET5M5+r- ( 4 - 5 ) 

1=0 i=l A i 

The index stands for the power of the first Chern class C\. The index j is 
the spin index, which takes the values j = 1, 2, ...,p — 1. 

We end up with an expansion, similar to the unitary case, but with dif- 
ferent coefficients 

i°gs=E<n^>n^ ( 4 - 6 ) 

For p = 3, the lowest orders of the O(N) model are given by 

iogz = ^E^i + t^(Ex2) 2 + 4^E^8 + 4^(E^i) 2 

+ s(E^)E^~C^ + o ( i) (4 . 7) 

Note that there is no odd-power of 1/A such as t - - This is due to the parity 
Aj — > — Aj for real antisymmetrix matrices of 0(2A r ). From the above series, 

the intersection numbers < Y\ nu j i T nij? > are obtained. In the unitary case, 
they are given by 

If s 

< r nun ■ ■■r naJa >= — _ n c i(A) n< C T (ii, • • -Js) (4-8) 

where c\ is the first Chern class and Ct{V) is the top Chern class [3]. In 
the present case, we call the intersection numbers as the coefficients of the 
expansion of logZ as (4.6). 

The numbers of r corresponds to the numbers of marked points s. The 
indices rij, ji are related to the genus g. 

E("i + -*-l) = (3-(l--))G/-l) (4-9) 
t=t P P 

The genus g is given through an expansion in powers of the inverse of the 
size of the matrix, as is standard for matrix models. For this purpose, we 
introduce an overall factor k in the exponent (3.5), an integral over k x k 
matrices. Then logZ/k may be expanded in a series in powers of k 2 ~ 29 , with 
genus g. For the present antisymmetric matrix Y, odd powers of are also 
present. (In the unitary case, only even powers appear). 
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The genus g is given by the Euler characteristics, 

V - E + F = 2 - 2(type) (4.10) 

where V, E and F are the numbers of vertices, edges and faces, respectively. 
For non-orientable surfaces, the genus g is replaced by the (type) : we can 
still use g but it takes half-integer values. This definition coincides with 
that of (4.9). In Fig.l-Fig.4, the lower order terms < r 10 > 9= i((l)torus 
and (2)Klein bottle), < Tq^ > g= i ((3) projective plane), and < r 2j i > g= i ( 
crosscapped torus) are depicted. 

For non-orientable surface , new characteristic terms are present, such as 
t\fi for the Klein bottle, in addition to the torus, and 1 (projective plane), 
t 2 ,i (crosscapped torus) which did not exist in the unitary case. From (4.7), 
we have the intersection numbers, 

1 2 1 1 2 1 

<r 1;0 > ff= i=-, < t 0i1 > ff= i= -, <r 2 , 1 > 3=| =— , <r lj0 > 9=0 =-. 

(4.11) 



Fig.l-Fig.4: opposite lines are glued together with matching arrow 




(1) torus 



(2)Klein 



(3)projective plane 




(4) crosscapped torus 



5 Evolution operators at edge singularities 



We have derived the higher Airy matrix model of (3.5) from the large N limit 
of the characteristic polynomials for the antisymmetric matrix Y, which is 
dual to the characteristic polynomial of X. The choice of the a n = a c = 
for (3.3) corresponds to a singular point in the spectrum of eigenvalues of 
X. At that critical point the density of states of X has a singularity at the 
origin. For a > a c there is a gap at the origin in the spectrum (whose support 
lies on the imaginary axis), and at a = a c this gap is closing. This happens 
also in the unitary case for an external matrix source with eigenvalues ±a 
[7, 8]. 

An integral representation for the evolution operators U for the vertices of 
O(N) (N even) may be obtained from the Fourier transform of the correlation 
functions. The derivations are given in appendix C. The evolution operators 
U(si, s n ) are defined as 

U( Sl , s n ) = ^ < tre SlX tre S2X • • -tre s " x > A (5.1) 

Fot the one point function of (2N) x (2N) antisymmetric matrix X, we have 
from (C.ll), 

U(«\ 1 *ir P ° x >- 1 l dV T\ V2 + (1 * ( V + l\ r sv + £ 

U{s) -2N <tTe >A --Wsf2ri}\(v + tf + a? U^ir 

(5-2) 

For a N x iV real antisymmetric matrix X, the sourceless probability 
density A = 0, 

JV(Af-l) 

P{X) = \e^\ Z={^) * (5-3) 
gives the expectation values 

< XijX k i >= ——(SikSji — SuSjk) 

< « 2 >= - N> \; 1) 

<(t,Y V >=^- 1 »< iV2 - iV + 4 ) 



I67 



2 



, N(N — 1)(2N — 1) , , 

< trf4 >= 1^ " ( 5 - 4 ) 

We have thus to compare 

U(s) = 1 < tre sX >= 1 + ^ < trX 2 > +^ < trX 4 > + • • • (5.5) 
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with this integral representation. For instance, in the case of N=l, (and 
7 = 1/2) the formula (5.2), after taking the residue at v — — |, leads to 



s 2 s 4 



U{ 8 )=e-r = l-- + - + ... (5.6) 



which indeed agrees with N = 2 in (5.4), 



l<trX»> = JJL^. = -i, 
N 2 ^=2 2 



1 <trA . a> = (N-1)(2N-1) 
N 4 



3 

N=2 4 



(5.7) 



For N = 2,3, it is easily verified that the integral representation of U(s) 
agrees with (5.4). 

In a previous article, we have found an explicit formula giving the zero- 
replica limit N — > for U(si, s n ) in the unitary case, in the absence of an 
external source (A =0) [5]. There it was shown that 

1 n s a 

hmf/( Sl ,..., Sn ) = -n2sh^- (5.8) 

n 

where a = s i- F° r n= l? this is simply 

i=i 

2 s 2 s 4 s s 

lim U(s) = 4§h^ = 1 + + -f-rr + • • • (5.9) 
n->o w s 2 2 24 5! • 2 4 v 1 

This means that 

lim — < trM 4 >= 1. (5.10) 
n^o N 

The replica limit counts the numbers of diagrams which can be drawn in 
one stroke line, and it corresponds to one marked point for the intersection 
numbers. 

For real antisymmetric matrices X, this replica limit of U(s) is different 
from the previous unitary case. From (5.2), we obtain 



dv , s . ( v + 



lim U(s) = -^-^-^(1 + ^)1— '^) e sv+ ^ 
n^o w s J 2m 2v \v + l) 

4 s 2 fT shx 
= — sh— + / dx 
Jo 



s 2 4 Jo x 

s 2 s 4 s 6 . . 

= 1 + — + — + + ••• 5.11 

4 96 1152 v ; 

The coefficients of s 2 and s 4 , \ and are consistent with (5.5) by (5.4), 
when we take 1 = \ and iV — > 0. The term of order s 2 is a Mobius band 



(projective plane), and it is a typical non-orientable surface. This term comes 
from the integral of (5.11), which does not exist in the unitary case. 

We obtain the connected part of the two-point correlation, after the shift 
Ui -> Ui + § in (C.19), 

lim U( Sl ,s 2 ) = - e ^ + ® l d ^l e ^^ 
X l °^ 1 + 2u 1 } [ui-(u 2 + fnul-(u 1 + fy] (5 - 12) 



with 



£/(si,s 2 ) = U(eis 1 ,e 2 s 2 ). (5.13) 

£l,e2=±l 



We deform the contour of u 2 near the origin to the poles (i) u 2 — u± + y, 
(ii)w 2 = —Mi — y, (iii)«2 = —Mi — 4f,(iv)w2 = Mi — y. Then we obtain in the 
replica limit N — > 0, 

hmC/(si, S2 ) = __L_ s h^( Sl + S2 )sh|( S i + S2 ) 

7 sh^-(si - s 2 )sh^(si - s 2 ) 



(si - s 2 ) 2 4 v zy 4 

+ \r +S2 d yU S -f^. (5.14) 
2 J Sl -s 2 1/4 4 

Since this is invariant under the change of signs of Sj, we have from (5.13), 

lim U(si, s 2 ) — 4 lim U(si, s 2 ). (5.15) 

AT->0 N^O 

For n > 2 we operate in a similar fashion and obtain the following result 
which generalizes the theorem for the unitary case [5] : 
Theorem 3: 



Urn U(s u ...,s„) = W(e 1 s 1 ,e 2 s 2 ,...,e n s n ), 

Ci=±l 

Sn ) = ^ ft(4sh^) + \ [ dy 1 - ft sh^ (5.16) 



1 A/.^^n , 1 r^lA.,^ 

where a = s x + ■ • • + s n . 



6 Intersection numbers from U(si, s n ) 

For one marked point, we consider the evolution operator U (s) in an external 
source A, chosen at a critical value. We discuss here the case p=3. From 
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(5.2), by the scalings of v — > yNv,s — > s/y/N, and by the critical value 
of = N, we have 

Exponentiating the term of power N, we have 

2 

-Mog[l + (v + ^) 2 ] + Mog(l + t» 2 ) + sv + ±j 

The first four terms are of order N after rescaling s — > TVs. Further, by the 
replacement by s — > y^s,^ — > 2iV,t> — ► u/y/2, we obtain in the large N limit, 

N 4 / \ 

U( ) = eTS f du Nsu^+^Ns^+Ns^u I U + S \ 

y ' Ns J 2m \u+ s ' 



2 . 



The shift u — > u — |, and w = -^t gives 



I /" dt Nt 3 + N s 8/3 t S 4/3 , 



We further make a scale t — * — it/(3N) 1 ^ 3 , then we have 

1 1 f°° t 3 
^) = 3V3 (iVg )W -y o ^os(-+xO 

/ dt-sm(- + xt)}. (6.5) 

7T 2 JO t 3 

where x = -iV 2/3 s 8/3 /(4 • 3 1/3 ). Using the Airy function Ai(x), 

\ f°° t 3 
Ai(x) = - / dtcos(—+xt) (6.6) 
n Jo 3 



we obtain 



1 / s 4 / 3 f3/VW 3 /■* 

U(s) = 3l/3(iVs)4/ 3 - l 2 J Jf ^'A(*')J (6.7) 

Since s is a Fourier transform variable, it is proportional to 

s ~ ^, * ~ s 8 / 3 ~ 1 (6.8) 
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The Airy function A^x) has asymptotic expansion 



Mx) = A(o)(i + + l -^x« + ■■■+ l -^p-x» + ■■■) 

+4(0)(* + + ^ + + ■ ■ ■) (6.9) 

where ,4,(0) = 3- 2 / 3 /r(2/3) and AJ(0) = -3- 1 / 3 /r(l/3). 

Therefore, we obtain the intersection numbers from the coefficients of the 
evolution operator U (s) for the one marked point like in the previous unitary 
case [5, 6]. The first series of expansion of A^x) in (6.9)gives the intersection 
numbers for the spin j — 1, and the second series gives the intersection 
numbers for the spin j = 0. From the first term in (6.7) with (6.9), one finds 
the intersection numbers < r n j > g for integer g 

i r(fi±!) 

< n8ff-5-i)/Sj > 9 = {24)9g] r( y ) (6-10) 

where j = for g — 31 + 1 and j = 1 for g = 3/. The intersection numbers for 
non orientable surfaces of half-integer genera, are obtained from the second 
term of (6.7), namely the integral of the Airy function. Taking into account 
the normalization (4.5), this leads to explicit results, such as 

1 

< r 2>1 > 9=3/2 - — 



7 Conclusion 

We have derived the intersection numbers for non-orientable surfaces from 
generalized Kontsevich Airy integrals over random antisymmetric matrices, 
the Lie algebra of the group SO(2N). An N — k duality bewteen k-point 
functions in N x N Gaussian matrix integrals, and N-point functions for 
k x ^integrals, in the presence of an external matrix source, allows one to re- 
late those generalized Airy integrals to the edge behavior of Gaussian models. 
Those Gaussian models are then much easier to deal with than the original 
integrals. The existence for Lie algebras of classical groups (such as Her- 
mitian matrices for U(N) or antisymmetric matrices for O(N)) of an Harish 
Chandra integral over the group elements is a key ingredient in these calcu- 
lations. These techniques should be useful for characterizing the geometric 
properties of those non orientable surfaces, and for comparing zeros of ana- 
lytic functions to random matrix spectra. 
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Appendix: derivation of theorem 1. 

Let us begin with the simple case, k=l. The determinant is given by the 
integral, 

< det(A • I - X) >=< J dc a dc a e Za{x - l - x)abCb > (A.l) 

over the 2N Grassmann variables c a ,c a . The probability measure P(X) for 
the average is 

p( X ) = ± e |trX2+trXA ^ 

Za 

Absorbing the antisymmetric part of the term c a X ab c b in the external source, 
A -> A', with 

A' ab = A ab - -{c a c b - c b c a ) (A. 3) 

We obtain 

trA' 2 = trA 2 - 2A ab c b c a - - (c a c a ) (c b c b ) (A.4) 

where a and b run over 1,2, ■•-,N. Writing the term (c a c a )(c b c b ) in the 
exponent as 

e -t^. i [ dye -^ (A . 5) 



/7T 

we obtain 

1 



< det(A • I - X) >= -j=jdyj dcdce- y2+{x+iy)SaCa - A « b ~ CbCa 

1 r N 

= -= dye- y2 T[det((\ + iy) ■ 1 - ^ia 2 ) 
V n J j= i 7 

1 r N 
= dye-y 2 H[(X + iy) 2 -a 2 ] (A.6) 

V 71 J j=l 



where 02 is the second Pauli matrix. On the other hand for the one-point 
function the right-hand side in (2.1) (theorem 1) is an integral over one single 
real number 



N 1 „ N , s 

<n<tot(« j .i-n> A = T Jn**{* 

1 r N 

= -^e- x2 / dyH(a 2 + y 2 )e-y 2 - 2iX y. (A.7) 



j= i 3=1 

N 
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The shift iy — > iy + A shows that (A. 6) and (A. 7) are identical. 

When k > 2 for averaging the k characteristic polynomials, k x 2N Grass- 
mann variables and c% are necessary (a — 1, k ; a, b — 1 • • • N). As for 
k — 1, we have 



trA — 2. !/„,'•„'•/, .''„''„'''//'/, :: : ,. ! ,, ' , ■ /■ 



1 

2 "a "a ""o ~o 2 



C„ C_ Ch Ch 



The last two terms are replaced by the following integrals, 



-£(e2eS)(<£<£) 



"^(^cf)(cfc«) = J dBe -^trB 2 +itrB^c^4 



(A.8) 
(A.9) 

(A.10) 
(A.ll) 



where B is a kx k Hermitian matrix and D is a kx k antisymmetric complex 
matrix. Thus we obtain, 



< ndet(A Q -X) >= f 

Q = l J 



xe 



\ jiffi I pa(3;(5 a I 1 n „/fi i 1 m 3 ' a 



(A.12) 



The exponent of this integral is a quadratic form in the Grassmann variables. 

Let us first consider the k = 2 case. The exponent is of the form 
En %M n V n , where 



- (r 1 c 1 r 2 r 2 r 1 r 1 r 2 r 2 ) 

^ n ~ \^2n+\i °2n+2> °2n+l5 °2n+25 °2ra+l> °2n+2' °2ra+l5 °2n+2/ 

with the 8 by 8 matrix M n 



where 



B r , 



D 



iB 12 




/ 


"^21 

V o 








-,021 

l -a n iB 2 i 
Ai + iSn 
A 2 + i-B 2 2 



D 21 \ 
D 21 





o y 





iB 21 
i . 
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12 



A 2 + iB 22 ) 



(A.13) 



(A.14) 



(A.15) 



(A.16) 
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Since the matrix M n is antisymmetric, the Gaussian integral over the \l/ n 
is the Paffian : 



Pf(M n ) 



\D 21 \ 2 - (Ai + iB u )(X 2 + iB 22 ) - \B 12 \ 
1 



7 



(Ax + iB n ) 2 - (A 2 + iB 22 ) 2 + 2|S 



12| 



2|£>2i| 
(A. 17) 



Writing 5 n = bi, B 22 = b 2 , B 12 = b 3 + 264, D — d\ + id 2 , we define the real 
antisymmetric matrix y as 



y 



61 64 + <i 2 &3 + di \ 

—b\ di — 63 64 — rf 2 

-64 - c?2 — di + h 6 2 

V — b 3 — di —64 + d 2 —b 2 



(A.18) 



J 



which satisfies the following equation, 



Since 



det(^ • I - Y) = Pf (M n ) 



try 2 = -2(b\ + b 2 2 )-A(bl + bl + d\ + dl) 
= -trB 2 -trD*D 



(A.19) 



(A.20) 



theorem 1 holds. 

For k > 2, the same procedure leads to the expression of an antisymmetric 
matrix Y, which is made of block of 2 x 2 matrices, given by 



Yij = (ImBij) ■ I + (ia^ReBij + ailmDij + <j 3 ReD 



'13 



(A.21) 



where the a,- are the 2x2 Pauli matrices. 



Appendix B: HarishChandra integral formula of theorem 2 

The Weyl group W of the SO(2N) Lie algebra is the permutation group 
S 2 n followed by reflection symmetries yi — > (e* = ±1), with an even 
number of ei = —1 . Then the sum over the elements of the Weyl group 
contained in the numerator of (4.1) becomes 

N 

1= E ( deta ) e M 2 J2 e -(j)y-u) x j] ( R1 ) 

e 1 =±l,---,e JV =±l ct£Sn j =1 
e 1 e 2 '"ejv = l 
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where the sum is restricted to reflections with an even number of sign changes, 
i.e. 

t l t 2 ---e N = l. (B.2) 

Since deter = (— (in which \a\ is the parity of the permutation), the sum 
over the N\ elements a of S N is a determinant. Therefore we obtain 

/= ]T det[e 2e ^'] (B.3) 

£1=±1,--',£]V = ± 1 
e 1 e 2 ---ejv = 1 

The result in this form is sufficient for the purposes of section 4. But one may 
go a bit further. Writing for each matrix element e 2€lViXj = cosh (tyiXj) + 
€j sinh (2yi\j) we obtain a sum of 2 N determinants weighted by products 
of 6j. The sum over those q, restricted by the condition (B.2), leads to a 
cancellation of all the terms except two . The final result is 

/ = 2 JV - 1 (det[cosh(2y i A i )] + det[sinh(2y i A J )]) (B.4) 

For instance, for N = 2 the signs are (ei, e 2 ) = (1, 1) or (—1, —1) . The sum 
over these terms gives 

/ = det[e 2yiX i] + det[e- 2yiX i] 

_ ^ e 2y 1 \ 1 +2y 2 \2 _ g2yi A 2 +2j/ 2 Ai ^ _|_ ^-2^1 Ai-2j/ 2 A 2 _ g -2yi A 2 -2j/ 2 Ai ^ ^ 

which is indeed identical to 

I = 2(det[cosh(2j/ i A j )] + det[sinh(2j/ i A j )]). (B.6) 



Let us apply the result (B.3) to the integral (3.5) : 

Z = J dYe-^ tlYP+1+tlYA (B.7) 

in which Y runs over the 2N x 2N antisymmetric matrices, the Lie algebra 
of SO(2N). One may use the rotational invariance of the measure to write 
Y = gyg^ 1 in which g is an element of SO(2N) and y is a canonical matrix 
(2.2), namely 



y = yiv © • • • © jjnv, v = ia 2 = 




The integral over Y may be replaced by an integral over g and over the y^s. 
The Jacobian is, up to a constant factor, J = Yl^jiyf — Uj) 2 - Using the 
Harish Chandra integral, one integrates over g and, using the result (B.3), 
one obtains 

Z= £ fdyi-- dy ^l " ^ e-^£^ +1 det[e 2 ^]. (B.8) 

e 1 e 2 ... £JV = l 
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Since p is odd, one can change — > yi. The antisymmetry of Y\{yf — yf) 
under permutations allows one to replace the determinant by its diagonal 
term exp(2J]yiAi) (up to a factor TV!). The restricted sum over the e; is 
simply the factor 2 N ~ l . We are thus led to the integral (4.2) of section 4. 

Appendix C: integral representations for U(si, • • • , s n ) 

For real antisymmetric matrices, the support of the density of states p(A) 
is the imaginary axis, 

p(A) = 1 < trS(X - X) > . (C.l) 

and it is an even function of A. The Fourier transform of p(A) is U{t), 

U(t) = J d\e ltx p{\) 

= < tre itx > . (C.2) 

Through an orthogonal transformation g (g G SO{2N)) ) one can bring the 
matrix X to the canonical form 

X = %\V © ■ • • © x N v, v — ia 2 — ^ j J ^ • (C.3) 

from which one has 

N 

tie sX = 2 cos(sXi) (C.4) 
i=i 

From theorem 2, the evolution operator U(s) becomes 

U(s) = ± E / nM^Tye- 2 ^^ 1 ' (C5) 

If we write cos(sa: a ) = Re(e lsXa ) the integral amounts to an integral over anti- 
symmetric matrices with the shifted source Oj — > dj = dj + \is5j ta . Then the 
integrals over the just the normalization for Gaussian antisymmetric 

matrices in the source A : 

J dXe L ^ x2+T * Ax = e-^ A \ (C.6) 
Indeed using again the Harish Chandra theorem, the left-hand side is simply 

J dXe^ x2+ ^ Ax = J X\d,- : ^0-, £*?+ 2 £»« (c.7) 
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which provides the result that we need : 

/ n^iA(x J V Ex?+2EwH " a = ^4re^ Tri2 ^ TrA2 . (C.8) 
This leads to 

1 N N ( (a 4- i^ 2 — ^j 2 \ 2 

= 4 £ n ( ( j' „. - + < s - < c - 9 > 

The normalization is such that U(0) = 1. It is useful to express (C.9) as a 
contour integral, 



7=1 \ U °7 / m 4 



where the contour encircles the poles u = a T 

Changing variables, u + y = —if, one obtains C/(s) 

^)-^/^n(^)^-*. (an, 

which is the representation that we have used in (5.2). 

The two-point correlation function U(s 1 , s 2 ) is given by Th.2, 

U( Sl ,s 2 ) = < tre SlX tre S2X > 



2 
N 



r A(x 2 ) 2 

J2 / ]Jdx i cos(s 1 x a Jcos(s 2 x a2 )-^-e~^ x ^ +2 ^ a ^. 

,a 2 =l i=l A v fl / 

(C.12) 



ai,a2 



We make replacements of cos(iSiX ai ) by \e lSiXa ^ (i = 1,2) in (C.12), and we 
name it as U(si, s 2 ). Then, we have 

U( Sl ,s 2 )= Yl U(e lSl ,e 2 s 2 ). (C.13) 

ei,e 2 =±l 

The sum in (C.12) is devided into two parts, J2 ai =a 2 an d J2 ai ^a 2 - The first 
part can be neglected. By the double contour integrals, the double sum is 
expressed as 

~ 1 r dudv " ( {u + ^) 2 -a 2 \ { (v + ^y-a^ 

U(Sl ' S2) = (Sri? { u ,_ a , ) { v ,_ a , ) 



(C.14) 



uv 

(iu - f)(iv - f ) [(u + l -f) 2 - v 2 ][u 2 -(y + ^f) 2 ] 6 
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By the Cauchy determinant identity, 

det^-^ = (-i)n(n-i)/2 ni<j(^-^)(y?-y?) (C15) 
A-v) UiAA-y]) 

with Xi = Ui + = Mj, above expression is simplified as 

xTTTT V 9 7 det 7 tt^ » (C16) 

For general n, we have similarly, 

N n /( , is±\2 _ 2\ i 

(M t + 2 S ^ J_ (C17) 



with 



C/(si, • • -,s n ) = ^(eisi, • • • ,e„s n ). (C.18) 

€4=±1 



This expression reduces to the previous one in (C.10) for n=l. We need the 
connected part of U(si, • ■ • , s n ), which is easily obtained from the expression 
of the determinant in (C.17). For n=2, we obtain by the change Ui — > —iiii, 

MM v «? + ; - f ) 2 - *2 - f ) 2 - • j 

This is the representation that we have used in (5.12). 
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